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ABSTRACT 


The Integrated Systems Program is capable of simulating both 
nominal and anomalous operation of the Apollo Cryogenics Storage 
System (CSS). Two versions of the program exist; one for the 
Apollo 14 configuration and the other for "J" Type Mission configur- 
ations. The program consists of two mathematical models which are 
dynamically coupled. A model of the CSS components and lines 
determines the oxygen and hydrogen flowrate from each storage tank 
given the tank pressures and temperatures, and the Electrical Power 
Subsystem (EPS) and Environmental Control Subsystem (ECS) flow 
demands. Temperatures and pressures throughout the components and 
lines are also determined. A model of the CSS tankage determines 
the pressure and temperatures in the tanks given the flowrate from 
each tank and the thermal environment. The model accounts for tank 
stretch and includes simplified oxygen tank heater and stratification 
routines. Supporting subroutines are used to calculate thermophysical 
properties from National Bureau of Standards (NBS) data and to calcu- 
late the CSS thermal environment during various mission phases. 
Anomalous conditions such as a fan and/or heater failure, degraded 
tank vacuum annulus performance, failed check values, and leaking 
lines can be simulated. The program is currently operational on the 
Univac 1108 computer. 
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INTRODUCTION 


The Integrated Systems Program is a mathematical representation 
of the Apollo Cryogenic Storage System (CSS) and is capable of simulating 
both nominal and anomalous operation of the system. The program is 
modular in structure and consists of a storage tank model, a compon- 
ents and lines model and auxiliary routines. Two versions of the 
program exist? one for the Apollo 14 spacecraft configuration and 
the other for the "J" Type Mission configuration. The "J" configur- 
ation differs from the Apollo 14 configuration primarily by the addition 
of a hydrogen storage tank and several other components. 

The tankage model is capable of simulating the thermodynamic 
performance of the CSS hydrogen and oxygen tanks. The model accounts 
for the effects of heat leakage, heater and/or fan operation, duty 
cycles, and control equipment as well as all significant thermal, 
elastic, thermoelastic and thermodynamic effects which influence tank 
performance. Capability to predict the effects of thermal stratifi- 
cation on the oxygen tank performance and to define oxygen tank heater 
temperatures are included. 

The components and lines model determines the oxygen and hydrogen 
flowrate from each storage tank based on the pressures and temperatures 
of the individual tanks and the flow demands of the Electrical Power 
Subsystem (EPS) and the Environmental Control Subsystem (ECS). A 
transient analysis is performed to determine the pressures and temper- 
atures at each junction and the flowrates in each leg of the system. 
These values are used to determine the instantaneous steady state 
pressure drop across all components in each leg. 

Anomalous operation of the tanks, such as heater and/or fan fail- 
ures, degraded tank vacuum annulus performance, and loss of tanks, may 
be simulated. In addition, line leaks and anomalous operation of 
any component may be simulated for the oxygen and/or hydrogen distri- 
bution systems. 

The auxiliary routines describe the thermal environment of the 
CSS and provide oxygen and hydrogen thermophysical data. The vehicle 
thermal model computes the outer wall temperatures of the storage tanks 
as well as spacecraft bay wall and shelf temperatures for various 
mission phases including translunar coast, transearth coast, attitude 
holds, and lunar orbit. The thermophysical data routine contains 
data for both oxygen and hydrogen obtained from the National Bureau 
of Standards (NBS) over the expected system operating range. 
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Input to the Integrated Systems program includes the initial state 
and quantities of the cryogens in the tanks, EPS and ECS flowrate demands, 
component characteristics, inertial accelerations and the mission mode. 

The program, written in Fortran V, is operational on the Univac 1108 
computer. 



SYSTEM DESCRIPTION 


The Apollo CSS supplies oxygen and hydrogen on a demand basis to 
the EPS and the ECS. The cryogenic fluids are stored in their super- 
critical states in tanks located in the Service Module (SM) . Each 
tank is equipped with individual components for filling; venting; 
quantity, temperature, and pressure measurements; pressure control 
and relief; and the accompanying displays and controls. For all 
missions through Apollo 13, the storage system was comprised of two 
hydrogen tanks and two oxygen tanks. The arrangement of these tanks 
in Sector IV of the SM is shown in Figure 1. For the Apollo 14 
Mission the storage, system consisted of an additional oxygen tank 
located in Sector I of the SM making a total of three oxygen tanks 
and two hydrogen tanks. Subsequent missions will have three hydrogen 
tanks as well as three oxygen tanks. The additional hydrogen tank 
will also be located in Sector I of the SM. Figure 2 presents the 
tankage arrangement in Sector I. 

During nominal operation, the storage system is maintained at a 
relatively constant pressure with fluid temperature increasing during 
the mission. Pressure is maintained above the critical pressure in 
the tanks by the application, either manually or automatically, of 
electrical heater power. Prior to the Apollo 14 Mission, electrical 
heat was supplied to the storage tanks in conjunction with destratifi- 
cation fans. However, subsequent to the Apollo 13 oxygen system 
failure, the system was redesigned to eliminate the fans in the oxygen 
tanks, although fans are still used in the hydrogen tanks. 

The nominal operating pressure of the storage tanks is 245 + 15 
psia and 900 + 35 psia for the hydrogen and oxygen tanks, respectively. 

The fluid flow rate is controlled on a demand basis by the ECS 
and EPS regardless of the storage system's pressure or operating mode. 
Under nominal operating conditions there are 29.3 pounds of hydrogen 
and 330.1 pounds of oxygen loaded into each of the storage tanks. 

Of these amounts, there are 28.2 pounds of usable hydrogen and 323.5 
pounds of usable oxygen. Other pertinent storage system operational 
parameters are presented in Table 1. 

Cryogenic Storage Tanks 

A major difference between the oxygen and hydrogen tanks is in 
the material used in their fabrication. These materials are listed in 
Table 1 and were obtained from Reference 1. In the broad sense, the 
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physical configuration for both tanks is similar in nature and they 
can be described as a typical unit. 

Each tank consists of two concentric spherical shells with 
thermal insulation in the annular space between the shells. In 
the oxygen tanks, the insulation supports the pressure vessel. The 
annulus of each tank is equipped with a vac-ion pump to maintain the 
low pressure required for adequate vacuum insulation, although in 
flight the oxygen pumps are the only pumps which may be operated. 

In addition, the pressure vessel supports (hydrogen tank only), fluid 
lines, and the electrical conduit are located within the evacuated 
annulus. The exit fluid line is routed through the annulus before 
exiting from the tank to provide vapor cooling. 

Prior to Apollo 13, each tank was equipped with a two element 
heater and two destratification fans. However, following the Apollo 
13 oxygen system failure, the two element heater and destratification 
fans in the oxygen tanks were removed and replaced with a three 
element heater. In addition, the heater thermostat for protection 
against high heater temperatures was removed and replaced with a 
temperature sensor. The new sensor is located on the heater tube. 
Another temperature sensor is located above the capacitance probe and 
measures the bulk fluid temperature. Electrical lines leading into 
the redesigned oxygen tank have also been sheathed to reduce the fire 
hazard. 

Cryogenic System Components and Lines 

Figures 3 and 4 are schematics of the component and line config- 
urations for the Apollo 14 oxygen and hydrogen systems, respectively, 
showing the relative locations of the components which make up the 
system. Figures 5 and 6 present similar schematics for the "J" Type 
Mission system. In general, the system Is composed of lines, valve 
modules, flow filters and flow restrictors. For missions prior to 
Apollo 14, the system valve modules for both the hydrogen and oxygen 
systems each contained two relief valves, two pressure transducers, 
two pressure switches and one check valve. For Apollo 14, a half 
valve module was added for oxygen tank 3 and contains one relief 
valve, one pressure switch, and a check valve. A solenoid operated 
shutoff valve was also added between oxygen tanks 2 and 3. 

The relief valves in these modules are a gradual opening 
differential pressure poppet type which assures the lowest possible 
leakage with virtually zero delta pressure between crack and full 
flow. These valves are designed to be unaffected by back pressure in 
the downstream plumbing and are temperature compensated over the full 
range of fluid temperatures which may occur. 
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The pressure transducers are absolute pressure devices which 
utilize bonded strain gauges to produce millivolt signals to the 
signal conditioner. The signal conditioner is integral to the 
transducer with a 0 to 5 VDC output linearly proportioned to tank 
pressure. The transducers can operate at environmental temperatures 
as low as -80°F and are capable of supplying a continuous readout of 
pressure over the specified range with a maximum error of + 2-1/2 
percent of full scale. 

The pressure switches are single-throw, single pole absolute 
pressure devices which may be used to control the motor switches 
which activate the destratification motors and/or tank heaters. These 
switches have been used successfully in an environment of -400°F. 

The check valve is a spring loaded poppet type and is designed 
to open at a differential pressure of approximately 1 psi. The valve 
seat has a large area to prevent chattering during flow in the normal 
direction and to help obtain a positive seal if pressurized in the 
reverse direction. 

The fuel cell valve module for the hydrogen system consists of 
two check valves and three solenoid operated valves in one valve 
body. The check valves are a dual seat type which allow the main 
passage to be full open and the auxiliary passage barely cracked at 
low flows. At high flows, both the main and auxiliary passages are 
completely open. The solenoid valves are poppet latch type and are 
actuated by a magnetic armature suspended on a Bellville spring. 
Opening and closing coils on both sides of the armature actuate the 
valve. The valves open against pressure and pressure helps to seal 
the valve against leakage when closed to flow in the normal flow 
direction. 

For the oxygen system, the fuel cell valve module contains check 
valves similar to those in the hydrogen fuel cell valve module. 
However, subsequent to Apollo 13, the solenoid valves have been taken 
out of the module and mounted externally and upstream of the flow 
transducers to the fuel cells. 

The oxygen and hydrogen inline filters consist of several 
chemically etched discs stacked in a cartridge. They are rated at 
5 microns nominally and 12 microns absolute with a capacity to contain 
.25 grams contaminant. 
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Several components have been added to the CSS for missions 
subsequent to the Apollo 14 ("J" type missions) . A half valve 
module containing a relief valve, a pressure switch, and a check 
valve has been added for the third hydrogen tank. There are three 
oxygen restrictors in the ECS for "J" type missions versus two for 
other Apollo missions. The third restrictor line is a branch of the 
oxygen tank 2 ECS supply line and reconnects downstream of the junction 
of the other two restrictor lines. This restrictor line includes a 
filter and eheck valve. 
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ANALYTICAL FORMULATION 


The Integrated Systems Program was developed as a tool to aid in 
the preflight planning and the postflight analysis of the Apollo Cryo- 
genic Storage System. Modular in structure, the program is comprised 
of a tankage model, a components and line model and several auxiliary 
subroutines. All significant effects which influence system perform- 
ance, such as thermal environment and mission mode, are included. 

Given the component characteristics, the initial state of the 
cryogens in the tanks, the EPS and ECS flow requirements and the 
mission mode, the program is capable of simulating the nominal perform- 
ance of the CSS. Figure 7 presents a flow diagram of the program. Out- 
put includes the temperature, pressure and flowrates at various points 
in the system as well as the remaining quantities and state of the 
cryogens in the tanks. 

In addition to nominal performance, the program is capable of 
simulating the operation of the CSS under anomalous conditions. The 
effects of degraded performance of the tank vacuum annulus caused by 
the presence of air or the cryogen in the jacket may be determined. 
Abnormal environmental temperatures or tank fan and heater failures 
may be simulated. Anomalous operation of the check valves, filters, 
relief valves or flow restrictors may be investigated. Oxygen or 
hydrogen line leaks may be modeled nearly anywhere in the system. 

Two versions of the program were developed. One version is 
capable of simulating the CSS up to and including the Apollo 14 config- 
uration. The other is capable of simulating the "J 11 Type Mission 
configuration. The primary difference in these programs is in the 
components and lines model. Only the "J" Type configuration will be 
discussed in this document. 

The following sections discuss the analytical formulation of the 
major models and supporting subroutines which comprise the Integrated 
Systems Program. 

Storage Tank Model 

The storage tank model was developed to simulate the thermodynamic 
performance of the Apollo CSS oxygen and hydrogen tanks. Given the 
initial conditions and flow demands, the model simultaneously updates 
the thermodynamic state in each storage tank based on control equipment 
operation, tank characteristics and heat leak. A flow diagram of the 
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model is presented in Figure 8. 

The model may be divided into two sections. The first section 
determines the maximum time increment over which all tanks may be 
simulated. The second section uses this time increment to update the 
state of the fluid in each tank corresponding with the conditions which 
prevail at the end of the time increment. 

The time increment depends on the rate of change of pressure and 
the operation of the pressure switches in each tank. The time for the 
tank pressure to reach a specified excursion (currently 15 psi) from 
its initial state is compared to the time required to reach the press- 
ure control limit. The shorter of these two time intervals is defined 
as the maximum dynamic simulation time increment for a given tank. 
Dynamic simulation time increments for all tanks are compared and the 
smallest one is chosen as the total tank system dynamic simulation time 
increment. This value is then compared with a user specified increment 
and the lesser is selected as the time increment for use in updating 
the state of all cryogenic tanks. 

Upon the determination of the time increment, the thermodynamic 
state of the fluid in each tank is updated. The average pressure for 
the time interval is determined from the initial tank pressure, the 
initial pressure rise (or decay) rate, and the length of the time inter- 
val. The average density for the time increment is determined from the 
tank volume, the initial weight of cryogen in the tank, and the flow- 
rate. Using these averaged thermodynamic properties, the average press- 
ure change rate is recomputed for each tank. The average pressure 
change rate and flowrate are then used to update the pressure and den- 
sity at the end of the time interval. Finally, fluid temperatures are 
obtained using the thermophysical routines. 

In addition to updating the state of each storage tank, the addi- 
tional flowrates required by the fuel cells to provide power for all 
electrical equipment operating during the time interval are determined 
and added to the EPS flow demand. The time required for a complete 
pressure cycle and heater cycle is also computed for each tank. 

One of the basic functions of the model is the determination of the 
rate of change of tank pressure with respect to time. The method of 
determining this derivative differs for the two cryogens since the 
effects of thermal stratification are considered significant in the 
oxygen tanks, but not in the hydrogen tanks. 

Oxygen Tank. The thermal stratification, resulting from non-uni- 
form heating of the cryogen, significantly affects the operation of the 
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oxygen tanks. The non-uniform heating leads to higher pressurization 
rates than that which would result from uniform heating. Furthermore, 
at high densities, thermally stratified oxygen can experience rapid 
pressure decay upon mixing due to the thermodynamic properties of the 
oxygen. 

In order to determine the pressure derivative, the oxygen tank is 
represented by three nodes or volumes. One node (V^) is associated 
with the fluid within the tank heater column. A second node (V 2 ) 
represents the fluid immediately surrounding the heater column. The 
final node (V3) represents the remaining fluid in the tank. The nodal 
representation of the tank is shown in Figure 9 . 

The oxygen which flows from the tank is withdrawn from V3. Mass 
interchange is permitted between V-^ and V 2 , but not between V 2 and V3. 

A result of this constraint is that V 2 continuously increases during a 
mission due to fluid withdrawal from V3. The effects of tank stretch 
and the compression of the fluid in the fill and vent lines (see 
Appendix II) on the pressure derivative are accounted for in V3. 

The cylindrical heater is described by a heat source of a specified 
temperature with a specified heat flux into the two adjacent oxygen 
nodes. The heater temperature data used in the model were obtained 
from Apollo 14 flight data and computer simulations of heater operation. 
The heat flux to the surrounding nodes is calculated from the electri- 
cal energy input to the heater, the heater configuration, and the time 
history of the heater temperature. Heat input is apportioned between 
the nodes adjacent to the heater by a simple conduction analogy which 
satisfies a total energy balance for the heater column. Heat conduct- 
ion is permitted between V2 and V3, but not directly between Vj and V2. 
The thermal energy exchange between Vj. and V2 must pass through the 
heater tube. All external heat leak is applied to V3. Heat leak radiant 
energy transfer is assumed to be negligible. 

Mass and energy balance relationships were developed for the fluid 
in each of the tank nodes. The energy balance equations are: 


d£ 

dt 

•Sl> 

II 

j q l + m l [ H + H ] 

U) 



-■1 <h r e i>i 


d£ 

dt 

= h j 
V 2 1 

| q 2 + m 2 j H ( m 2^ h l + H ( -m 2^ h 2 ] 
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" p 2 ®2 " WV 


344 



and 


d£ _ 43 

dt v 3 


dV 3 • 

q 3' p 3 0 3 dT~ + m 3 



where: 


H(x) 


0 x < 0 

1 x > 0 


(3) 


dV3 = dV t ank d£ _ ^ _ dV L 

dt dp dt dt dt 

The mass balance equations are: 



m 3 - "withdrawal 

The Nomenclature used in these equations is defined in Appendix I. 

The mass balance relations may be substituted into Equations (1) , 
(2) and (3) to yield a set of three non-linear equations in the three 

variables m, and -^2 which may be solved iteratively. The result- 

ing pressure derivative is then applied as described above. 

The pressure collapse potential is defined as the difference 
between the instantaneous value of pressure and the value that would 
result from sudden complete mixing of the fluid. Using the individual 
values of enthalpy for the three volumes , a mean enthalpy is determined 
for each tank. The mean density of the fluid in each tank is then 
determined from the known weight of the remaining oxygen and tank 
volume. The tank pressure for this mixed condition is obtained from 
the thermophysical routines. Finally, the pressure collapse potential 
is computed using this pressure and the pressure under stratified 
conditions . 

Hydrogen Tank . During the Apollo 7 Mission, flight tests were 
performed which demonstrated that thermal stratification did not 
significantly affect the performance of the hydrogen tanks. It was 
therefore assumed that thermodynamic equilibrium existed in the hydro- 
gen tanks at all times. The equation which determines the rate of 
change of pressure with respect to time for the hydrogen tanks is: 
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(4) 


Equation (4) is derived in Appendix II and accounts for variations in 
tank volume due to fluid pressure and temperature excursions as well 
as the compression of fluid in the fill and vent lines. All terms of 
Equation (4) may be determined from the CSS system and fluid thermo- 
dynamic properties. 


Components and Lines Model 

The components and lines model simulates the operation of the 
Apollo CSS from the tanks to the user subsystems. Given the EPS and 
ECS flow demands, and the tank conditions, the pressures, tempera- 
tures and flowrates are determined at various points in the system. 

A flow diagram of the model is shown in Figure 10. 


The model is divided into two parts which are called, for 
convenience, the matrix and nodal models. The matrix model is a 
mathematical representation of the oxygen and hydrogen systems based on 
the linearized electrical analogs shown in Figures 11 and 12, respectively. 
A transient, compressible flow analysis of each system is performed to 
determine the flow rate in each leg and the pressure at each matrix 
node. The nodal model then uses this information to perform a steady- 
state, incompressible flow analysis of each leg to determine the 
pressure at selected locations along that leg and to update the non- 
linear characteristics of the individual components and line segment 
resistances. The nodal model also performs a transient thermal analysis 
of the lines and components. Detailed electrical analogs of the oxygen 
and hydrogen systems are shown in Figures 13 and 14, respectively. 


Matrix Model . The conservation of mass principle is used to 
determine the matrix node pressures. For any node i, this may be 
written (see Appendix III) 


El . 


l 


dE. 


dt 



(5) 


346 



where 


and 


V 

i 


**±. 1*1 





( 6 ) 


( 7 ) 


The matrix model solves this set of simultaneous equations for each 
system by an implicit backward differencing technique. The flow rate 
in each leg is then evaluated using the conductances in the character- 
istic matrix. 


The program can simulate up to two leaks in each system at any 
location except downstream of relief valves. Adding a leak is equi- 
valent to adding a matrix node. These are included in the character- 
istic matrix as floating columns and rows. The flow at leaks is 
assumed to be choked and is evaluated in the same manner as over- 
board vents which are discussed below. 


Nodal Model . The nodal model performs a steady state analysis 
of each leg using the matrix model node pressures and temperatures as 
boundary conditions. The individual components and line segments are 
analyzed as follows : 


Lines and Fittings . The pressure drop in a line segment 
from inlet 1 to outlet j is calculated by (see Appendix III): 



Since fluid velocities are small, it is assumed in the derivation of 
Equation (8) that the flow is laminar and incompressible. The con- 
servation of energy principle is used to determine the line tempera- 
tures. It is assumed that conduction through the stand-off brackets 
and along the lines is negligible. Then, for any node, k. 
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( 10 ) 


, , dT 

'"S’k dT 


°V 


< T S * 


T k >- UA 4T LM 


Fittings have an equivalent L/D of approximately 3. Bends are 
included by an equivalent L/D which is a function of the bend 
angle (Figure 15) . Clogged lines are simulated by incorporating 
an increased equivalent length into Equation (8) . 

Filters, Check Valves, Relief Valves, and Shutoff Valves . 

The pressure drops across filters and all open valves are evaluated 
from Equations (8) and (9). An equivalent L/D of 100 is assumed 
for filters, open check valves, and open relief valves. The equi- 
valent L/D for open shutoff valves is assumed to be 10. The 
pressure drop across all closed valves is calculated from: 



where the conductance G, is assumed to be 10 ^ for check valves 

and 10 ^ for relief valves. Check valve and relief valve 
hysteresis is modeled as shown in Figures 16 and 17, respectively. 
Clogged filters and valves are simulated by incorporating an equi- 
valent L/D that is larger than the nominal value. Failed check 
valves and relief valves are simulated by incorporating a cracking 
and/or reseat pressure that will hold the valve open or closed as 
desired. 


The temperature of filters and individual valves are deter- 
mined from 

dT 


<m Vk tt ■ - UA AT LM, 


( 12 ) 


while the temperatures of valves contained in a valve module 
housing are determined from 


(me ) 


dT k 


p k dt 


V A 

= (^L ) (T,-T,.) - UA AT 


1 k 


LM 


(13) 


Restrictors . The pressure drop across the oxygen system flow 
restrictors is evaluated using the data in Figure 18. It is 
assumed that the process is isothermal. 


Overboard Vents . For compressible flow such as overboard 
venting, the flow is choked. The pressure is evaluated by: 
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(14) 



P i A.K 


The nodal model also checks the solution for convergence. The 
updated resistance of each component and line segment is compared to 
the resistance used to calculate the flow rates. If the difference 
in resistor values is within the convergence criteria, the program 
proceeds to the next time step; if not, the flow rates are re-evalu- 
ated based on the updated resistances. Iteration of the characteris- 
tic matrix proceeds in this manner until the solution converges (nor- 
mally within two iterations). 


Auxiliary Routines 

In order to adequately simulate the operation of the Cryogenic 
Storage System, it was necessary to develop several auxiliary routines. 
Two of the most important routines will be discussed in this section. 

They are the vehicle thermal model and the thermophysical properties 
routine. 

Vehicle Thermal Model . The vehicle thermal model was developed 
to compute the Cryogenic Storage System environmental temperatures and 
to determine the heat leak of the oxygen and hydrogen tanks under 
normal conditions or under conditions of loss of vacuum in the insula- 
tion annulus. The heat leak into the storage tanks is primarily a 
function of the flowrate through the vapor cooled shroud (VCS) and the 
temperature in the bay. The bay temperature depends to a large extent 
upon whether the vehicle is in an inertial hold mode or a passive thermal 
control mode. The effect of vacuum loss in the jacket is to increase 
the heat leak into the tank. 

Oxygen Tank Heat Leak . The heat leak into the oxygen tanks 
during nominal operating conditions is computed in the thermal 
model using empirical data correlations developed by the Beech 
Aircraft Company (Reference 2). These correlations relate the 
heat transfer into the oxygen tank with the temperature of the 
oxygen in storage, the oxygen expulsion rate and the tank outer 
wall temperature. These correlations are given by the following 
equations: 

Q = AL1 (T 2; ~* - T 2 * 5 ) - m (h - h) (15) 

r 1 vcs vcs K ' 
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(16) 


Q = AL2 (T 2-5 - T 2 ' 5 ) 
r vcs 2 

q = AL3 (T x - T 2 ) + .9Q r 

The value of the constants which apply to Equations (15), (16), 
and (17) are as follows: 

Redesigned Tank: ALl = .9205 x 10 _5 BTU/HR°R 2,5 
AL2 = .7278 x 10 3 BTU/HR°R 2 * ^ 

AL3 = .03775 BTU/HR°R 


The amount and temperature of the oxygen in storage is 
computed in the tankage model of the Integrated Systems Program 
and is input to the thermal subroutine. 

The outer wall tank temperature is computed within the ther- 
mal subroutine through a series of data curve fit correlations. 

The data for these correlations were obtained using the Thermal 
Mathematical Model (TMM) of the Apollo Service Module. This 
program is composed of over 800 lumped parameters that constitute 
a thermal network specifically formated for solution with the 
SINDA Multi-Option computer program. Using TMM, a set of para- 
metric computer runs was performed varying the mission mode 
(inertial hold, passive thermal control or lunar orbit), the 
quantity of the cryogen in the tanks, and the flowrates from the 
tanks. A uniform internal tank temperature was assumed in deter- 
mining the temperature of each oxygen and hydrogen tank surface. 

For use in the thermal subroutine, the temperatures computed 
for the surface of the oxygen tanks by the TMM were averaged for 
each tank and curve fit as a function of time. The form of the 
resulting curve fit equations is: 

T = a Q + a^x + a 2 x 2 + a 3 x 3 (18) 

where x is the time argument. The coefficients for the polynomial 
curve fit equations are shown in Table 2. Using these equations, 
the tank outer wall temperatures are first estimated as a function 
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of mission mode, time in mode and quantity. These values are then 
adjusted for the sun look angle. The temperature of the fluid 
in the VCS is then computed from the temperature of the fluid in 
storage and the outer wall temperature using an iteration scheme. 
Once the temperature of the VCS fluid has been determined, the 
heat flux is computed using Equations (15), (16), and (17). 

Hydrogen Tank Heat Leak . The heat leak to the hydrogen tank 
under nominal operating conditions is computed in the thermal 
model using data taken from the Apollo Flight Support Handbook 
(Reference 1). This data, presented in Figure 19, is a correl- 
ation of estimated heat flux as a function of the fluid expulsion 
rate and the quantity of fluid remaining in the tank. This data, 
however, does not include the effect of temperature in the bay on 
the heat flux. Since the outer wall temperature is not taken into 
account, the mission mode cannot be considered in the calculation 
for the hydrogen tank heat leak. In addition, the discontinuity 
occurring in the curves at a flow rate of 0.25 lbm/hr does not 
appear reasonable. Although the correlations presented in Figure 
19 do not appear to adequately define the heat leak into the 
hydrogen tanks, it is the only data currently available. However, 
if additional data becomes available, the hydrogen tank heat leak 
correlations used in this portion of the program may be updated 
easily. 

Vacuum Annulus Failure . The nominal pressure in the tank _ ( 
annulus under the design vacuum condition is approximately 2 x 10 

psia (10 ^ mm Hg) . Depending upon the type of failure, the 
annulus pressure in the failed condition will be somewhere between 
the nominal pressure and the annulus burst disc pressure. Burst 
disc pressures for the oxygen and hydrogen tanks are approximately 
75 and 80 psia, respectively. Within these pressure bounds, gas 
in the annulus may exist in either rarefied or continuum regimes. 
The criteria for determining the gas regime is based primarily on 
the ratio of the molecular mean free path to the amount of space 
available for molecular motion. This ratio is defined as the 
Knudsen number (K n ) . For values of Kj < 0.01 the gas is normally 
treated as continuum flow, while for values of > 1.0 the gas 
is assumed to be in the rarefied regime. For the CSS tanks, an 

_ O 

annulus pressure on the order of 10 J psia represents the trans- 
ition between the two regimes. 

The type of analysis used to compute the amount of heat 
transferred across the annulus depends upon the flow regime. For 
the rarefied regime the equation for determining the heat flux 
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per unit area is given by Reference 3 as: 


q = A„ 


(— ) 
'n-1' 


£. 

8n v/mt" 


(T 2" T l ) 


(19) 


For continuum flow regimes , the heat flux per unit area is given 
by Reference 4 as : 


q = -kr 


(T^) 

(r 2 -r 1 ) 


( 20 ) 


An analysis was performed for the Apollo cryogenic storage 
tanks (Reference 5) assuming residual gas in the annulus consisted 
of air or the appropriate cryogen. The analysis showed that for 
an annulus pressure of less than approximately 10~7 psia, heat 
conduction by the residual gas was essentially negligible. However, 
heat conduction increased rapidly with increasing annulus pressure. 
The total heat transfer to the tanks as a function of the annulus 
pressure is shown in graphical form in Figures 20 and 21. These 
curves have been incorporated into the thermal model such that 
when the thermal subroutine encounters a failed tank flag, the 
data presented in the figures is used to determine the heat flux 
to the tank based upon the vacuum jacket pressure. This pressure 
is an input to the program. 


Components and Line Environmental Temperatures . The environ- 
mental temperatures used to determine the state of the fluid in 
the components and lines were obtained in a manner identical to 
the tank environmental temperatures. The Thermal Mathematical 
Model of the Apollo Service Module was used to obtain data at 
several locations as given in Table 3. These data were incorpor- 
ated into the Integrated Systems Program as a series of curve fit 
correlations. The form of the resulting equation is the same as 
Equation (18). The coefficients for the curve fit equations are 
presented in Table 2. 


The TMM did not model the Command Module (CM) cabin, the CM 
Reaction Control System Tank Annulus, or the ECS umbilical. The 
temperature at these locations was assumed to be 70°F. 


Thermophysical Properties for Oxygen and Hydrogen . Computer sub- 
routines for computing the thermophysical properties for oxygen and 
hydrogen were developed and included into the Integrated Systems 
Program. The majority of the data used in these subroutines was 
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recommended by the National Bureau of Standards (NBS) . That data not 
available from the NBS was taken from Reference 6 (Stewart). In general 
these subroutines determine pressure (p) , temperature (T), density (p), 
enthalpy (h), the specific heat at constant pressure (c ), the specific 
heat at constant volume (c ) , viscosity (y) , thermal conductivity (k) , 
isotherm and isochor derivatives and the functions <|> and 0, providing 
sufficient information concerning the fluid state is known. The data 
is included in the subroutines both in equation and table form. A 
routine is also included in the subroutines for interpolating between 
values from the tables. 

Oxygen Properties . Data contained in the subroutines for 
computing the thermophysical properties of oxygen were taken from 

References 6, 7, and 8. Reference 6 presents an equation of state 
for oxygen determined by R. B. Stewart. The isotherm and isochor 
derivatives were derived from this equation. While Reference 7 
(Weber) also presents data for the isotherm and isochor derivatives, 
the results derived from Reference 6 are valid for a greater range. 
Both sets of data were included in the thermophysical subroutines 
to provide greater flexibility since the derivatives derived from 
Reference 6 may be determined with known values of temperature 
and density, while the data from Reference 7 uses either pressure 
and temperature or pressure and density. 

The equation of state provided by Reference 6 may be used to 
compute pressure from known values of density and temperature. 

This equation is valid for temperatures between 117 and 540°R and 
for pressures to 5000 psia. The equation of state is as follows: 


M K J. 


+ (n 6 T 2 + n ? T + n g + n 9 /T + n^/T^p 3 ’ 

+ (n ll T + n i2^ 4 + (n i 3 + n i4 T) e 5 
+ p 3 (n 15 /T 2 + n 16 /T 3 + n^/T 4 ) exp (n^p 2 ) 

r n _ 


+ 


+ P 5 (n 18 /T 2 + n 19 /T 3 + n^/T 4 ) exp (n^p 2 ) 

+ P ? (n 21 /T 2 + n 22 /T 3 + n^/T 4 ) exp (n^p 2 ) 

W^V 28 -*; 28 ) exp [^ e Xn) 2 ] (21) 
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where : 


P c - 13.333 gmol/liter, critical density 
T c = 154.77 °K, critical temperature 


In this expression, p is in atmospheres, T is in °K, and p is in 
gmol/liter. The constants for this equation are given in Table 4. 
The subroutine in which Equation (21) is programmed converts the 
units of density from lbm/cu-ft to gmol/liter and temperature 
from °R to °K prior to use. The output pressure is in psia. 


The isotherm derivative, 


lE. 

3p I T, 


is the partial derivative of 


Equation (21) with respect to the density at constant temperature 
and is given by: 


3p 


= R T + 


24 

E 

i=l 


n . X . 
l i 


( 22 ) 


Values for the variables in this equation are presented in Table 
5. 


The coefficients n^ are the same as those for Equation (21). The 
subroutine which computes the isotherm derivative converts density 
from lbm/cu-ft to gmol/liter and temperature from °R to °K prior 
to its use in Equation (22). The output derivative is in units 
of psia-f t 'Vlbm. 

3 1 

The isochor derivative is the partial derivative of 

O T | P y 

Equation (21) with respect to temperature at constant density and 
is given by: 

24 

= PR + £ n -, x i (23) 

i=l 

where the coefficients n^ are the same as those listed for Equation 

(21). The variables in this equation are shown in Table 6. The 
conversion from one set of units to the other is performed the same 
as for Equations (21) and (22) . The output derivative is in urtits 
of psia/°R. 


l£ 

3T 
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In addition to providing data for the isotherm derivative and 
the isochor derivative, Reference 7 presents data for density, 
enthalpy, specific heats (at constant pressure and volume) and the 
functions 0 and <j>. All are in tabular form. Further flexibility 
was provided by permitting these properties to be determined from 
known values of pressure and temperature or from pressure and 
density. These tables are valid for temperatures from 100 to 540°R 
at 5 degree increments and for pressures ranging from 367 to 1250 
psia at increments of 73.5 psi. Table 7 lists the thermophysical 
parameters which are available. 

Data from Reference 8 is also in the form of tables and 
provides viscosity and thermal conductivity as functions of 
pressure and temperature. This data is valid for temperatures 
ranging from 180 to 540° R at increments of 18 degrees and for 
pressures from 367 to 1250 psia at increments of 73.5 psi. 

Hydrogen Properties . The thermophysical properties used 
in the hydrogen subroutine are based upon computer programs 
developed by the NBS (Reference 9). These properties are valid 
for pressures from 1 to 5000 psia and for temperatures ranging 
from 25 to 5000 °R. In the NBS programs all thermophysical 
properties are determined from known values of pressure and temp- 
erature and include enthalpy, density, specific heats for both 
constant pressure and constant volume, viscosity and thermal con- 
ductivity. However, in the subroutine developed for the Integra- 
ted System Program, provision was made to compute the pressure and 
temperature from other known properties as described below. The 
pressure, density and temperature are in units of psia, lbm/ft , 
and °R, respectively. 

The NBS program expresses hydrogen density as a polynominal 
equation in terms of pressure and temperature, p = p (p,T). The 
equation has been programmed such that pressure may be determined 
from this equation as a function of density and temperature using 
a Newton-Raphson iteration method. An a priori value of pressure 
is required to initiate the iteration and should be estimated as 
closely as possible to the actual value. Temperature is computed 
in a manner similar to that of pressure from this same expression. 

The available parameters are listed in Table 7. 
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TABLE 1 


CRYOGENIC SYSTEM OPERATIONAL PARAMETERS 



Hydrogen 

Oxygen 

TANK MATERIAL 

5 Al-2.5 Sn ELI tl 

Inconel 718 

TANK WEIGHT (PER TANK) 

. 


Empty (Approx) 

80.00 lb. 

90.82 lb. 

Usable Fluid 

28.15 lb. 

323.45 lb. 

Stored Fluid (100% 
indication) 

29.31 lb. 

330.1 lb. 

Residual 

4% 

2% 

Maximum Fill Quantity 

30.03 lb. 

337.9 lb. 

TANK VOLUME (PER TANK) 

6.80 ft 3 

4.75 ft 3 

TANK PRESSURIZATION 



Heaters elements per tank 2 

3 

Flight Resistance 

Nominal Voltage 
Power 

78.4 ohms per 
element 
28 V DC 
10 watts per 
element* 

18.45 ohms per 
element 
28 V DC 

42.5 watts per 
element* 

Total Heater 
Heat Input Per tank 

68.2 BTU/Hr 

434.8 BTU/Hr 

Ground Resistance 

78.4 ohms per 
element 

18.45 ohms per 
element 

Voltage 

Power 

28 V DC 
10 watts per 
element* 

28 V DC 
42.5 watts per 
element* 

Total Heater 
Heat Input Per Tank 

68.2 BTU/Hr 

434.8 BTU/Hr 


357 




TABLE 1 (Continued) 



Hydrogen 

Oxygen 

Pressure Switch 



Open Pressure Max. 

260 psia 

935 psia 

Close Pressure Min. 

225 psia 

865 psia 

Deadband Min . 

10 psia 

30 psia 

Destratification Motors 
(2 Motors Per Tank) 



Voltage 

115/200 V 

None 


400 cps 


Power - Average 

3.5 watts per 
motor* 

None 

Total Average 
Motor 

Heat Input Per Tank 

23.8 BTU/Hr 

None 

SYSTEM PRESSURES 



Normal Operating 

245 + 15 psia 

900 + 35 psia 

Spec Min. Dead Band 

10 psi 

30 psi 

of Pressure Switches 



Relief Valve Note: Relief Valves are referenced to environmental 

pressure, therefore pressure at sea level 

(psig) will be the same 

value In vacuum (psia) 

Crack Min. 

273 psig 

938 psig 

Full Flow Max. 

285 psig 

1010 psig 

Reseat Min. 

268 psig 

865 psig 

Outer Tank Shell 
Burst Disc 

+10 


Nominal Burst Pressure 

90 -20 psid 

75 + 7.5 psid 

SYSTEM TEMPERATURES 



Stored Fluid 

-425 to 80 °F 

-300 to 80 °F 

*Conversion Factor: 1 watt = 3.41 BTU/Hr 
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TABLE 2 


CURVE FIT COEFFICIENTS REPRESENTATIVE 
OF THE CSS THERMAL ENVIRONMENT 


Mission 

Mode 

Node 

a 

0 

CF) 

a l 

(°F/Hr ) 

32 2 
(°F/Hr ) 

a 3 

(°F/Hr 3 ) 


Time A 

rgument 

T 

inf 

°F 

t+3 

(Hrs) 

■nai 

SH 

PTC 

Tanks 1 & 2 

46. 152 

35. 715 

105.97 

-38.106 




PTC 

Tank 3 

25.164 

73. 761 

202.92 

-112.08 



X 

Hot Soak 

Tanks 1 & 3 

60.0 

4.44 

.0 

0.0 


X 


Hot Soak 

Tank 2 

58.523 

31.715 

-3.57 

0.1277 


X 


Cold Soak 

Tanks 1 & 3 

-15.824 

368.40 

-485.5 

232.26 



X 

Cold Soak 

Tank 2 

-61.812 

549. 39 

-650.84 

283.34 



X 

PTC 

150 

66.0 

0.0 

0.0 

0.0 

66.0 


X 

PTC 

151 

85.3 

-64.0 

0.0 

0.0 

80.0 


X 

PTC 

152 

47.0 

0.0 

0.0 

0.0 

47.0 


X 

PTC 

153 

49 .5 

28.6 

0.0 

0.0 

57.0 


X 

PTC 

154 

42.5 

46.5 

0.0 

0.0 

48.0 


X 

PTC 

155 

117.3 

-130.0 

0.0 

0.0 

101.0 


X 

PTC 

156 

132.4 

-145. 3 

0.0 

0.0 

118.0 


X 

PTC 

157 

91.8 

-77.3 

0.0 

0.0 

86.0 


X 

PTC 

158 

70.0 

0.0 

0.0 

0.0 

70.0 


X 

FTC 

159 

48.4 

18.2 

0.0 

0.0 

49.0 


X 

PTC 

163 

47.0 

0.0 

0.0 

0.0 

47.0 


X 

Cold Soak 

150 

80.0 

-3.9 

0.27 

0.0 

56.0 


X 











TABLE 2 (Cont'd) 

CURVE FIT COEFFICIENTS REPRESENTATIVE 


OF THE CSS THERMAL ENVIRONMENT 


Mission 

Mode 

Node 

a 

o 

(°F) 

a i 

(° F/Hr ) 

a 2 

(°F/Hr 2 ) 

a 3 

(°F/Hr 3 ) 

T 

inf 

°F 

Time 

Arcument 

t+3 

(Hrs) 

t+3 

(Hrs) -1 

Cold Soak 

151 

115.2 

-18.6 

1.4 

-0.036 

31.0 


X 

Cold Soak 

152 

109.3 

-22.3 

1.7 

-0.041 

5.0 


X 

Cold Soak 

153 

88. 7 

-10.2 

0.45 

-0.0074 

7.0 


X 

Cold Soak 

154 

107.5 

-18.3 

1.2 

-0.03 

6.0 


X 

Cold Soak' 

155 

60.7 

-0.25 

0.046 

-0.0021 

58.0 


X 

Cold Soak 

156 

51.3 

7.65 

-0.42 

0.0067 

96.0 


X 

Cold Soak 

157 

99.6 

-12.6 

1.07 

-0.029 

50.0 


X 

Cold Soak 

158 

50.0 

5.8 

-0.26 

0.0039 

96.0 


X 

Cold Soak 

159 

111.6 

-20.1 

1.42 

-0.035 

5.0 


X 

Cold Soak 

163 

32.1 

12.2 

-0.70 

-0.013 

109.0 


X 

Hot Soak 

150 

103.3 

-10.2 

0.812 

-0.002 

50.0 


X 

Hot Soak 

151 

47.7 

8.73 

-0.443 

0.00071 

103.0 


X 

Hot Soak 

152 

32.1 

12.2 

-0.7 

0.013 

109.0 


X 

Hot Soak 

153 

105.2 

75.8 

-882.4 

396.6 

105.0 


X 

Hot Soak 

154 

97.2 

62.1 

-870.3 

395.0 

99.0 


X 

Hot Soak 

155 

120.7 

40.4 

-975.4 

450.0 

122.0 


X 

Hot Soak 

156 

24.6 

17.9 

-0.892 

0.013 

136.0 


X 

Hot Soak 

157 

wm 

17.6 

-1.2 

0.028 

105.0 


X 



CURVE 


Mission 

Mode 

Node 

Vi 

a 0 

(°F) 

Hot Soak 

158 

86.1 

Hot Soak 

159 

102.2 

Hot Soak 

163 

109.3 

Lunar Orbit 

150 

53.5 

l.unar Orbit 

151 

74.6 

Lunar Orbit 

152 

45.0 

Lunar Orbit 

- 153 

79.0 

Lunar Orbit 

154 

58.0 

l.unar Orbit 

155 

51.1 

Lunar Orbit 

156 

58.9 

Lunar Orbit 

157 

69.0 

Lunar Orbit 

158 

76.72 

Lunar Orbit 

159 

57.0 


Lunar Orbit 


163 


45.0 
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a + a.X J - a_ X : + a~X' where X = time argument, 
o I Z J 





TABLE 3 


ENVIRONMENTAL TEMPERATURE LOCATIONS 


SINK INDEX 

LOCATION 

150 

TUNNEL 

151 

Bay 4 

152 

Bay 4 

153 

Forward Bulkhead 

154 

Bay 4 

155 

Bay 4 

156 

Bay 4 

157 

Bay 4 

158 

Forward Bulkhead 

159 

Forward Bulkhead 

160 

Umbilical 

161 

CM RCS Tank Annulus 

162 

CM Cabin 

163 

Bay 1 


362 




TABLE 4 


CONSTANTS FOR STEWART'S EQUATION OF STATE 


R = 

0.0820535 



n l5 = 

1.73655508 

X 

10 




-3 




5 

n l = 

3.38759078 

X 

10 

n 16 = 

3.01752841 

X 

10 

n 2 = 

-1.31606223 



n 17 = 

-3.49528517 

X 

10 7 

n 3 = 

1.92049067 

X 

10 3 

n 18 = 

8.86724004 

X 

io _1 

n 4 = 

1.92049067 

X 

10 7 

n 19 = 

-2.67817667 

X 

10 2 

n 5 = 

-2.90260005 

X 

O 

t-H 

0 

n 20 = 

1.05670904 

X 

10 5 




-8 




-3 

n 6 = 

-5.70101162 

X 

10 0 

n 21 = 

5.63771075 

X 

10 

n 7 = 

7.96822375 

X 

io “ 5 

II 

CN 

CM 

C 

-1.12012813 



n 8 

6.O7O22502 

X 

10“ 3 

n 23~ 

1.46829491 

X 

10 2 








-4 

n 9 = 

-2.71019658 



n 24 = 

9.98868924 

X 

10 

n l0 = 

-3.59419602 

X 

10 

n 25 = 

-0.00560 



n ir 

1.02209557 

X 

10- 6 

n 26 = 

-0.157 



n 12 = 

1.90454505 

X 

10-* 

n 27 = 

-0.350 



n 13 = 

1.21708394 

X 

10- 5 

n 28 = 

0.90 






-3 





n i/. = 

2.44255945 

X 

10 
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X 1 

= 2pT 

X 2 

= 2p 

X 3 

= 2p/T 2 

X 4 

= 2p/T 4 

X 5 

= 2p/T 6 

„ 2 2 

X 6 

= 3p T 

X 7 

= 3p 2 T 

X 8 

= 3p 2 

X 9 

= 3p 2 /T 

X 10 

= 3p 2 /T 

X 11 

= 4p 3 T 

X 12 

= 4p 3 


= f„f f- 

24 

9 2 3 

exp 

(n 25 P 2 ) 

n 28 

n 28 

P 

" P c : 

ex P [ n 26 f 2 + 

2n 27 

XT-I c )£ 3 

2f x p 

“25 


X 13- 5(> 

x 14 = 5p 4 /t 

X 15 

X 16‘ f 6 /?3 
X 17 V T ‘ 
X 18 * f 7 /l2 

X 19 V" 3 

X 20 " f 7^ T 

X 21- V 2 
X 22 - f 8 /T3 
X 23 * V 1 ' 


<n 28+l ) n (, '28+l ) f 

f 10 f 3 + P f 2 f ll 


2 3 

= 3f jp £ + f 5 p J 

= 5f lP 4 + f 5 p 5 
= 7f lP 6 + f 5 p ? 
^ n 28+l^ p 28 

= n 2g p (n 28-i> 


f ll 2f 2 n 26 f 3 f 10 




TABLE 6 


VARIABLES FOR COMPUTING ISOCHOR DERIVATIVE 
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TABLE 7 


THERMOPHYSICAL PROPERTIES SUBROUTINES 


PROPERTY 

UNITS 

INPUT 

VARIABLES 

REFERENCE 

OXYGEN 

Density 

lbrn/ft"^ 

P&T 

7 

Enthalpy 

btu/lbm 

P&T 

7 

Enthalpy 

btu/lbm 

P&p 

7 

Function 0 

btu/lbm 

P&T 

7 

Function <j> 

3 

psia-ft /btu 

P&T 

7 

Isochor 

psia/°R 

T&p 

6 

Isochor 

psia/°R 

P&T 

7 

Isochor 

psia/°R 

P&p 

7 

Isotherm 

3 

psia-ft /btu 

T&p 

6 

Isotherm 

3 

psia-ft /btu 

P&T 

7 

Isotherm 

3 

psia-ft /btu 

P&p 

7 

Pressure 

Psia 

T&p 

6 

Specific Heat, 
C P 

btu/lbm-°R 

P&T 

7 

Specific Heat, 
C P 

btu/lbm-°R 

P&p 

7 

Specific Heat, 
C 

V 

btu/lbm-°R 

P&T 

7 
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TABLE 7 (CONT’D) 

THERMO PHYSICAL PROPERTIES SUBROUTINES 


PROPERTY 

UNITS 

INPUT 

VARIABLES 

OXYGEN 

Specific Heat, 
C 

V 

btu/lbm-°R 

P&p 

Thermal Conduc- 
tivity 

btu(f t/hr) / 
ft 2 /°R 

P&T 

Temperature 

°R 

P&p 

Viscosity 

lbf-hr/ft 2 

P&T 


REFERENCE 


HYDROGEN 



Density 

lbm/ft^ 

P&T 

9 

Enthalpy 

btu/lbm 

P&T 

9 

Function 0 

btu/lbm 

P&T 

9 

Function 4> 

3 

psia-ft /btu 

P&T 

9 

Isochor 

psia/°R 

P&T 

9 

Isotherm 

3 

psia-ft /btu 

P&T 

9 

Pressure 

Psia 

T&p 

9 

Specific Heat, 
C 

P 

btu/lbm- °R 

P&T 

9 

Specific Heat, 
C 

V 

btu/lbm-°R 

P&T 

9 

Thermal Conduc- 

btu(ft/hr)/ 

P&T 

9 

tivity 

ft 2 /°R 



Temperature 

°R 

P&p 

9 

Viscosity 

lbf-hr/ ft 2 

P&T 

9 
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FIGURE 1 . SM SECTOR IV CSS TANKAGE ARRANGEMENT 
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FIGURE 2. SM SECTOR I CSS TANKAGE ARRANGEMENT 
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FIGURE 3 

,L0 14 CRYOGENIC STORAGE SYSTEM 
OXYGEN SCHEMATIC 
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FIGURE 6 

APOLLO J TYPE CRYOGENIC STORAGE SYSTEM 
HYDROGEN SCHEMATIC 
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FIGURE 7 


INTEGRATED SYSTEMS PROGRAM 
FLOW DIAGRAM 
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FIGURE 8 


STORAGE TANK MODEL 
FLOW DIAGRAM 
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FIGURE 9 

OXYGEN TANK NODAL MODEL 
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FIGURE 10 

J TYPE COMPONENT AND LINES 
MODEL FLOW DIAGRAM 
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FIGURE 13 

ELECTRICAL ANALOG OF THE OXYGEN SYSTEM FOR THE NODAL MODEL 
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F - Filter 








































FIGURE 18 

RESTRICTOR PRESSURE DROP* 
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HYDROGEN TANK HEAT LEAK 
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APPENDIX I 
NOMENCLATURE 

Area , in^ . 

Accommodation coefficient, equal to 1.0. 

V 4 2 

Capacitance analog, C = * lbm - in /lbf 

ZKip, 

Specific heat of fluid at constant pressure, btu/ 
(lbm-°R) . 

Specific heat of fluid at constant volume, btu/ 
(lbm - °R). 

Wall thickness, in. 

Line diameter, in. 

2 2 4 

Electrical potential analog, E = p , lbf /in . 
Energy of system 

Friction factor, dimensionless. 

2 4 

Conductance analog lbf -hr/ (lbm- in ) 

2 

Gravity constant lbm-in/ (lbf-hr ) 


Enthalpy of stored fluid, btu/lbm 

Enthalpy at initial state, btu/lbm 

Enthalpy of fluid in vapor cooled shroud, btu/lbm. 

Current analog, I = m, lbm/hr. 

Residual current analog, I = ~ lbm/hr. 

Thermal conductivity, btu/ (hr-f t-°R) . 

Constant , 


LUUUULL1V 1 L_y 9 UtU/ U XX/ • 

:, K =(~~) ' " 2(5=1), lbm-°R 1/2 / (lbf-hr). 


Knudsen number, A/1, dimensionless. 
Characteristic length, ft. 
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L 

M 

m 

m 

n 

N 


N 


RE 


P,P 


q 

q 


Length, in. 

Gas molecular weight, lbm/ (lb-mole) 

Mass , lbm. 

Mass flowrate, lbm/hr. 

Ratio of specific heats, Cp/c^, dimensionless. 

Poly tropic exponent for gas lines attached to tank, 
dimensionless. 

P V D 

Reynolds number, — - — , dimensionless. 

Pressure, psia 

Total heat input rate, btu/hr. 

2 

Heat leak flux, btu/ (hr- ft ) 



r 


R 



s 


t 


T 

T 

T 

T 


1 

2 


vcs 


u 


Radiation heat leak, btu/hr. 

Mean pressure vessel radius, in. 

Universal gas constant. 

Inner tank wall radius, ft. 

Outer tank wall radius, ft. 

Fluid velocity, in/hr. 

Time, hr. 

Temperature, °R. 

Temperature of tank outer wall, °R. 

Temperature of tank inner wall, °R. 

Average temperature of fluid in vapor cooled shroud, °R. 
Internal energy of fluid, btu/lbm. 


390 



u 

v 

V 

V 1 

V 

a 

(3 

A 

e 

0 

X 

V 
P 
a 

v 

♦ 

* 


2 

Overall heat transfer coefficient, btu/ (hr- in -°R). 

3 

Specific volume, ft /lbm. 

3 

Tank Volume, ft . 

3 

Volume of lines attached to tank, ft . 

Young's modulus for tank material, psi. 

Tank material thermal expansion coefficient, 1/°R. 
Tank thickness, in. 

Difference, dimensionless 


Emissivity, dimensionless 

Thermodynamic function, -P -jr | , btu/lbm. 

P 

Mean free path, ft. 

Viscosity, lbm/ (in-hr). 


3 

Fluid density, Ibm/ft . 

Stephan-Boltzman constant, = 1.1x10 ^ , btu/ (hr-in 
°R 4 ). 


Poisson's ratio for tank material, dimensionless. 

JL 2ei 3 

Thermodynamic function, p 8u 1 psia-ft /btu. 
Line bend angle, degrees. 


391 



Subscripts 


c 

i 

j 

k 

1 

LM 

r 

R 

s 


Cracking 

Inlet 

Outlet 

Arbitrary Line or node 

Arbitrary line or node 

Logarithmic mean 

Reseat 

Residual 

Sink 
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APPENDIX II 


DERIVATION OF dp/dt FOR VARIABLE VOLUME TANKS 
DUE TO PRESSURE AND TEMPERATURE CHANGES 


Introduction . 

An expression for the rate of change of pressure for the cryo- 
genic storage tanks has been developed for a constant volume 
(Reference A-l) and for a variable volume caused by tank stretch due to 
pressure changes 1 (Reference A-2). The latter expression is: 

(AII-1) 


-u) 


(AII-2) 
(AII-3) 

If Young's Modulus, Y, is taken as infinity (no tank stretch), Equation 
(AII-1) reduces to the expression for dp/dt as presented in Reference 
A-l. 


d£ 

dt 


*0 ‘ , i 
v ■ m + v 

1 + 30^>p.r 
1 2dY 


I 


where : 


0 = - 


. 3h 
I ( iE, 

p '3u ; p 


The expression for dp/dt considering variable volume resulting 
from tank stretch due to both pressure and temperature changes is 
developed in this Appendix. For completeness, a portion of the 
derivation from Reference A-2 is also presented. 

Derivation . 

The First Law of Thermodynamics may be written: 


dEn 

dt 


q + m h' 


dV 

P dF 


If it is assumed that the system is at uniform pressure, p, 
and that kinetic and potential energy may be neglected, then 


En = mu = pVu 
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and 


Now 


d( pVu) 

dt 


q + m h' - 


dV 

} dt 


ud(pV) + p V^ = q + m h' - p d ffi^ 
dt 


Then 


u = h - pv. 

udm + pVd (h-pv) _ + ^ h * _ P^dv _ pvdm 


dt 


dt 


dt 


dt 


urn + p V dh - pV d(pv) = q + m h'- pm dv - pvm 
dt dt dt 

° • 

m (u + pv) + pV dh - pVp dv - pVv dp = q + m h' 
dt dt dt 


m h + pV dh - m y/av - V dp = q + m h' - p yf&v 
dt y' dt dt / 

• 

V dp =p V dh - q + m (h - h') 
dt dt 

• 

dp = p dh - Q_ + m 
dt dt V V (h - h'). 

Since hsh', the last term is negligible. 


Thus , dp = p dh - cj_ 
dt dt V 


Assume that 


h = h(p,p) 


then 


dh = 9h dp + 3h 
3p p 9p 


dp 


dv 

pm IF 


(AII-4) 
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dh 

dt 


3h 

3p 


dp + _3h d£ 


p dt 


3p 


dt 


Now, 


and 



dp 

dt 


.(£) _ 1 dm 

V- " V dt 

dt 


m ' dV £ dm £ dV 

v 2 dt m dt " V dt . 


So , 


dh = 

3h 1 d£ 

+ 3h 

i 

p dV ] 

dt 

3p 1 p dt 

3p 

1 p 1_ m dt 

V dt J 


Substituting for© from Equation (AII-2) 



dh 

dt 

= Jh 
3p 

d£ - 
p dt 

Since 

h = 

u + 

pv = u 

then 

3h 


in 1 + 


3p 

P 

3p 1 P 


3h 


in I + 


3p 

P 

3p 1 P 

So , 

dh 

_r 

isil + 


dt 


3p 1 P 

Substituting 

for 0 from 


dh 


'I + 


dt 


P<p 


0 m + 0 dV 
m V dt 

+ *> 

P , 



1 

P 


ll dp 0 ’ 0 dV 

p J dt m m V dt . 

Equation (AII-3) 

1] dp 0 * 0 dV 

p J dt m m V dt . 



(AII-5) 
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It is at this point that the assumption is made that tank stretch is a 
function of both pressure, p, and temperature, T. 


dV _ _3V I d£ _3V dT 

dt 3p 1 1 dt 3T p dt . (AII-6) 


It is necessary to determine dT/dt 
dT _ ill _3T I _dp 

dt 3plp dt 3p Ip dt . (AII-7) 

Now, m = pV 

dm dp tT , dV 
dt = dt V +p dt 


So , d£ m £ dV 

dt V " V dt (AII-8) 

Substituting Equations (AII-8) into Equation (AII-7), and substituting 
Equation (AII-7) into Equation (AII-6) results in: 

dV = 3V I d£ + 3VI rjTI d£ + 3Tj (| - £ ^-)1 

dt 3p It dt 3T|p |_3p| P dt 3plp V V dt J 


For convenience, the following substitutions have been made: 


A = — 

A 3p It 


B = 3V 
B 3T Ip 


C -|1 , 

3p I P 


Then, 


I 

3p I p 


dV = A d£ L d£ f H _ dV \ 1 

dt A dt dt ° ( V V dt'J 


fr (A + BC)^ + Br,S-Saf 
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dV 

dt 


(1 + 


BDP . 

v ; 


(A + BC) 


dp , BDm 
dt V 


dV 

dt 


(A + BC) 


1 + 


dp BDm 

dt V 
BDP 
V 


(AII-9) 


Substituting Equation (AII-9) 


dh Jl ll d£ 
dt [ P< P p -1 dt 


into Equation (AII-5). 

e r < A + BC > $ * 60 7 ] 

+ V L , , BDP J 

V 


• • 

dh _ [ 1_ 1 Q(A + BC) I dj3 0m ©BDm 

dt L P<f> P V+BDp J dt " m V (V+BDp) . 


But from Equation (AII-4) 
dh _ _1 d£ . jL_ 

dt p dt pV 


Therefore, 

y * • 

Iffi + l. JL , l/ 0 (A+BC) 1 d£ 0 m 0 BDm 

c/ dt pV [ P4> + /T + V+BDp J dt m V(V+BDp) 

and , 

r 1 0 ( A+BC ) I dp _ 0m 0BDm + £_ 

^ p V+BDp J dt m V(V+BDp) pV 

resulting in: 


d£ 

dt 


Qm 

PV 


1 

p? 


. q _ 0 BDm 
PV V (V+BDP) 


a .jA+JC) 

V+BDp 
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or, 

• • 

On $ J>£_ QdiBDmo 

d£ = V V V (V+BDp ) 

dt 60o (A+BC) 

V+BDp 


Substituting for A, B, C and D 

0<J>mP ^ L 


Qm6 

V 


6q 

v 


TT 


3T 

p 3p 


iE 

dt 


, 3v 
v(v+ t> W 


3T 

P ^P 


1 + 


06p(. 


3V 

<>E_ 


+ iv 

T 3T 


3T 

£.i£- 


) 


„ . 3V 
V + p'rrr 
w 3T 


_3T 
P 3p 


(AII-10) 


Reference A-2 presents the derivation for change in volume with respect 
to pressure at constant temperature as: 


3 P IT (1-u) . (AII-11) 


The change in volume with respect to temperature may be obtained as 
follows : 


3V 

3T p 


3V 

3r 


_3r 
p 3T 


P 


3V 

3r 


3_ 

3r 


( I m 3) " 


4 irr 


2 _ 3V 
r 


For small changes in temperature, then 


3V 3V Ar 

3T p r AT 


3aV 


(AII-12) 
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Substituting these expressions into Equation (AII-10) results in: 


Qm4> 

V 

„ _ 

dt 


+ 13. _ 

V 


' . 3T 

0<J>m 3ap -t— 

3p P 

3T I 

V (l+p3a3p |p) 


1 +9 * p <Hy (1 - HSL > 


1 + 3ap 

3p p 


(AII-13) 


It is necessary to account for one additional effect which results in 
a variation in the change in volume with respect to pressure. This 
effect is the compression of the fluid in the fill and vent lines of 
the tank. These lines extend from the pressure vessel to the warm 
environment. The fluid at the cold ends is relatively incompressible 
as compared to the near ideal gas conditions at the other ends. When 
pressure increases, the fluid in the lines is compressed by the 
expanding stored fluid which move up the lines. To account for this 
effect, it was assumed that the lines will be filled with an ideal gas 
at ambient temperature and that no mixing between the tanked cryogen 
and the perfect gas would take place. Then: 


pV = constant. 

JLj 


Thus , 


U n -dp _l_ dV T 

V L dt + npV L ~ 
dt 


= 0 


or 


The mass flow, 


dV 

dt 

m, 


L _ dp 

np dt 

is the total flow from the tank: 


(AII-14) 


m = m + m_ (AII-15) 

o 2 

where: m^ = flow to EPS and ECS subsystems, lbm/hr 

m = p L = flow from fill and vent lines, lbm/hr 

1 dt 
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Note: Flow out of the tank is negative in sign. 

Substituting Equation (AII-14) into Equation (AII-15) yields: 

0 V L 

m = m - p — -rr 
o np dt 

When this equation is substituted into equation (AI 1-13) and rearranged, 
the resulting expression is: 


0 £ 

V 


d> 

m o + 


0 £ 

V 


' , 3T 

m 3ap-r— 
o 3p 

(1 + p 3a 


3T , 
3p P J 


dt 


1 + 


0£ pV Q± 

V " V 

np 


V L 3ap 


3T 

3p 


np (1 + 3ap 


3T 

3p p 


+ F 

) 


(All-16) 


where: 


F 


,3r 


3T 


Q4>P ^2dY (1 ~ U) + 3a 3^ 

1 + 3ap lr 

3p p 



Thus, Equation (AII-16) describes the rate of change of pressure for 
a variable tank volume considering both pressure and temperature 
changes. Note that if the coefficient of thermal expansion, a, is 
equal to zero and if the plumbing fluid compression volume, V L , is 
zero. Equation (AII-1) is obtained. Table A-l presents the 
constants used for both the oxygen and hydrogen tanks. The partial 
derivatives are obtained by differentiating the equation of state. 
Since the Integrated Systems Program uses an equation of state that 
is a function of temperature and density, it is desirable to obtain 
3T/3p in the following form: 


3T 

3£ 

JR 

3p 

P 3P 

T 3T 


This relation was taken from Reference A-4, page 32, Equation (3 11). 
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Then 


3T 

3p 


1 


i£ 


3p 

T 3T p 


_ iRi ill 

3p 1 T 9p 1 p 


Thus , this partial derivative may also be obtained from the equation 
of state. 
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TABLE A-l 

CRYOGENIC TANK PROPERTIES 


OXYGEN 

Material Inconel 718 

3 

Volume, V, ft 4.75 

Inside Radius, r, in 12.513 

Young's Modulus, Y, psi 30 x 10^ 

Poisson's Ratio, u, dimensionless 0.29 

Wall Thickness, d, in 0.059 

Coefficient of Thermal Expansion, a, /°R 5.93 x 10 ^ 


HYDROGEN 

Material 

Volume, V, ft^ 

Inside Radius, r, in 

Young's Modulus, Y, psi 

Poisson's Ratio, u, dimensionless 

Wall Thickness, d, in 

Coefficient of Thermal Expansion, a, /°R 


Titanium- 5AL-2 . 5SN ELI 
6.80 
14.103 
17 x 10 6 
0.30 
0.044 

1.23 x 10" 6 


Note: Values for a were derived from Reference A-3, Table 7.3.4, 

page 7-59. 
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APPENDIX III 


DERIVATION OF EQUATIONS USED IN THE 
COMPONENT AND LINE MODEL 

MATRIX MODEL 

The matrix model determines the pressures for each of the oxygen and 
hydrogen system nodes shown in Figures 11 and 12, respectively. 

The conservation of mass principle may be written: 

I m = v f£ • (AIII-1) 

For a perfect gas, 


p = pRT , 


(AIII-2) 


Combining Equations (AIII-1) and (AIII-2) and performing the 
differentiation results in: 


l m = -X- (A - pgi). 

RT 2 dC ^ 

To set up an electrical analogy, let 

E = p2 

for an equivalent driving function, and 

I = m 

for the current. 


(AIII-3) 


(AIII-4) 


(AIII-5) 


In terms of the transforms variables, Equation (AIII-3) becomes: 


I 1 = 


V dE v/j dT 

2RT» / E dt RT 2 dt 


or 


1 1 - # - 

L dt 


where 


C = 


V 


2RT/E 


(AIII-6) 


(AIII-7) 

(AIII-8) 
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(AIII-9) 


is an equivalent capacitance, and 

= dT 

R RT 2 dt 

is an equivalent residual current. 


NODAL MODEL 

The nodal model determines the pressures of the oxygen and hvdroeen 
along each leg of the oxygen and hydrogen, Figures 13 and 14, 
respectively. 


The instantaneous pressure drop for incompressible flow through a 
line with inlet i and outlet j is given by the Darsy-Weisbach equation. 


/4fL N ps^ 

p i - ■ < T' ) . . W 

J 11 C 

For incompressible flow in a circular tube, 

7TD 2 pS 

m = pAs = t . 


(AI.II-10) 


(AIII-11) 


Combining Equations (AI1I-10) and (AIII-11) and simplifying results in: 

8m 2 


P i P j D ) . j g c r 2 D 4 p ’ 


(AI1I-12) 


Let the mean density between i and j be based on the average pressure 
and temperature. Then, 


ij 


1 p.+p. 

= — i l 
R — — 


T.+T. 
i 1 


(AIII-13) 


Substituting Equation (AIII-13) into Equation (AIII-12) and simplifying 
yields : 


2 2 

P- ~ P. = 
1 1 


8R(T +T. ) 
1 1 

g 7T 2 D 4 
B c 


(ML) 

V D ’ 


ij 


nr 


(AIII-14) 


For laminar, incompressible flow in a circular tube, 

16 4iryD 

1 " N DU ' • . (AIII-15) 

RE m 

Combining Equations (AIII-14) and (AIII-15) and simplifying results in: 
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Pl 2 - pA Wl * 


(AIII-16) 


In terms of the electrical analog variables, Equation (AIII-16) becomes: 


E i - E J ' G 


where 


— 2 - 8 ^ -- (T +T . ) 


(AIII-17) 


(AIII-18) 


An instantaneous heat balance through a line with inlet i and out- 
let j is: 


where 


UA AT t „ = m c (T.-T.) . 
LM p 3 i 


( vv - <vv 


V T i 


(AIII-19) 


T -T 

k j 


(AIII-20) 


Combining Equations (AIII-19) and (AIII-20) and simplifying yields: 


T = 1 - e p T. + 


( 4c Jl . , "(*0 


P T . (AIII-21) 
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